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We report on the calculation by the MILC Collaboration of the electromagnetic effects on kaon 
and pion masses. These masses are computed in QCD with dynamical (asqtad staggered) quarks 
plus quenched photons at three lattice spacings varying from 0.12 to 0.06 fm. The masses are fit 
to staggered chiral perturbation theory with NLO electromagnetic terms, as well as analytic terms 
at higher order. We extrapolate the results to physical light-quark masses and to the continuum 
limit. At the current stage of the analysis, most, but not all, of the systematic errors have been 
estimated. The main goal is the comparison of kaon electromagnetic splittings to those of the 
pion, i.e., an evaluation of the corrections to "Dashen's theorem." This in turn will allow us to 
significantly reduce the systematic errors in our determination of m„/m,(. 
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Introduction. - The disentangling of electromagnetic (EM) and isospin- violating effects in the 
kaon and pion systems is a long-standing problem. Understanding these effects is crucial for com- 
puting light quark masses, which are fundamental parameters in the Standard Model and important 
for phenomenology. Indeed, the size of the EM contributions to the kaon masses is the largest un- 
certainty in determinations of m u /m c i from the lattice [p, and in particular in our calculations [g]. 
The contributions have until recently been taken from a variety of phenomenological estimates, and 
therefore have quite large and not well controlled errors. We have been working on reducing these 
uncertainties for some time by calculating the EM effects directly on the lattice; progress has been 
reported previously in Refs. [^J |J. 

The error in m u /m c i is dominated by the error in the mass difference {M\± — M^y, where y 
denotes the total EM contribution, i.e. , the difference between the value of a quantity in the presence 
of electromagnetism, and its value in a world in which all EM charges, both of valence and of sea 
quarks, are turned off. One may try to relate this difference to the much better understood difference 
in the pion system, (M^ ± — M^ o y. To lowest order (LO) in chiral perturbation theory (ZPT) these 
EM splittings are in fact the same; this observation is known as Dashen's theorem [|0. We aim to 
calculate on the lattice the corrections to Dashen's theorem, which may be parameterized by 

(M 2 K± -JKfo) 7 = (1 + e)(M^ ± -M'- o y . (1) 
Our computations employ full QCD but quenched photons. As pointed out in Ref. [^], however, the 
EM-quenching effects on e may be calculated and corrected to NLO in XPT, with controlled errors. 
Similarly, squared mass differences of the form (Mp —Mp,y are calculable with controlled errors in 
our setup, where P is any light pseudoscalar meson and P' is the corresponding meson made from 
neutral valence quarks with the same masses as those in P. We also compute the EM effects on the 
K° alone, namely (M^ o y. In this case, however, the quenching effects are not calculable in XVI, 
and uncontrolled errors remain. For m u /nid, the uncertainty coming from (M^. ) 7 is, fortunately, 
subdominant. 

In the pion system, isospin-violating effects on the mass splitting are known to be small 
(see, e.g., Ref. [|I|]), so the experimental splitting is almost completely electromagnetic: (M^ ± — 
M^,) expt f» (M^± —M^y. It would be costly to simulate the true 7T°, which has EM disconnected 
diagrams even in the isospin limit. Instead, we simulate a "7To" whose squared mass is a simple 
average of the squared masses of uu and dd mesons, computed with connected diagrams only. Be- 
cause all EM contributions to neutral mesons vanish in the chiral limit, both the true (M^ o y and our 
(M^^.y are small in any case. From Zweig-rule considerations, we suspect that the disconnected 
contribution is smaller still. We estimate 

(M^y~(M^,-Ml,y. (2) 

Since there are errors coming both from quenched electromagnetism and from the neglect of dis- 
connected diagrams, one might expect ~ 100% errors in this estimate, however. On the other hand, 
our calculation of the EM pion splitting, using (Ml ± -M^y « (M\ ± -M^ „y suffers only from 
the neglect of disconnected diagrams, not from uncontrolled quenching effects. As a rough estimate 
of the former we take 50% of (M^y, calculated through Eq. @. 

Chiral Perturbation Theory. - We fit our lattice data to expressions from (partially quenched) 
staggered chiral perturbation theory (SZPT) in order to extrapolate to the physical light quark 
masses and to the continuum. We consider Goldstone (taste pseudoscalar mesons composed of 
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valence quark x and valence antiquark y, with masses m x and m y . Let M xy $ be the mass of such 
a meson with valence-quark charges q x and q y (q xy = q x — q y is the meson charge), and let AM^, 5 
be the squared-mass splitting AM^, 5 = M 2 5 — M 2 Y , 5 , where the primes in the second subscript 
indicate that the valence-quark charges are set to zero. To NLO in SZPT, AM 2 5 is given by 

<2 I * 2 



AM^ 5 = q 2 xy 8 E M-^ 1 e 2 q 2 xy Jt(l 5 ^31n(^, 5 /A^) - 
2,3 1 

"167^2 16 ^ [* Ew9 *V^.« ta (^,S ) " lyalxr^yo,!; H-^yo,S, \ 
+c\q 2 xy a 2 + c 2 <^ (2m/ + m. v ) + C3 {q 2 x + (m^ + m y ) + c^q 2 ^ {m x + m y ) + c 5 (^m* + q 2 m y 



(3) 



where 5sm is a low-energy constant (LEC), a runs over the sea quarks, t, runs over the staggered 
tastes, Cj are the LECs at NLO, and m/, ra s are the light and strange sea-quark masses. At this order, 
the meson masses denoted by ^# on the right hand side may be taken to be the tree-level masses 
in the absence of electromagnetism. Note that Ref. [^] mentions an additional NLO analytic term 
proportional to {q\ + q 2 )a 2 ; this is not possible for a Goldstone splitting because it does not vanish 
as m x + m y for small quark masses in the limit q y = q x , where EM effects are chirally symmetric. 

Our statistical errors in AM^, 5 are ~0.3% for charged mesons and ~ 1.0% for neutral mesons. 
It is clear that NLO SZPT cannot be expected to give a good description of the splittings at that 
level of precision. For reasonable fits, NNLO terms are needed. The SZPT logarithms have not 
been calculated at that order, but we add all possible analytic terms. This may be justified by noting 
that the NNLO logarithms will be small at low mass, where the extrapolation is performed, and will 
be well approximated by analytic terms in the region near m s . 

To allow for finite-volume (FV) effects, we include standard terms dependent on m n L s from 
EM tadpoles (L s is the spatial lattice size), as well as an empirical EM finite-volume correction 
of the form f v q 2 y /L 2 used previously in Ref. [Q], where f v is a constant. However, as discussed 
below, our measured finite- volume effects are rather small at present, and including or omitting the 
finite-volume terms from the fits makes little difference in the final results. 

Lattice setup. - We calculate the meson spectrum in quenched EM backgrounds on a set of 
asqtad ensembles with 2+1 flavors and 0.12 fm>a>0.06 fm. See Ref. [Q] for a table of lattice 
parameters. The valence quarks have charges ±2/3e, ±l/3e or 0, where e = e p h ys , 2e p h ys , or 3e p h ys 
(e p hy S is the physical electron charge). For the results reported here, only the e = e p h ys data is used. 

Figure shows some of our partially quenched data for the splittings of physically charged 
and K mesons. We investigate the FV effects on two ensembles with a 0. 12 fm, ami = 0.01, and 
am s = 0.05, but different volumes (L = L s /a = 20 and 28). The figure shows a comparison of the 
FV effect seen in our data with the difference between the two volumes expected from the results of 
the BMW Collaboration [§]. The effect we see is smaller: 0.35(45) of the difference expected from 
BMW. However, our larger-volume ensemble has only 274 configurations, leading to the rather 
large errors in the comparison. We are currently increasing the statistics on that ensemble in hopes 
of clarifying the issue. 

We have measured the taste splitting between the Goldstone pion and the other local pion on 
our ensembles. The amount of taste violation caused by high-momentum photons is nonnegligible 
for mesons made from quarks with higher-than-physical charges, although it is always significantly 
smaller than that caused by high-momentum gluons. The fit function Eq. (|3]> is based on the neglect 
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of taste violations caused by photons; that is the reason that we focus here only on the data with 
physical quark charges. Given that photon-induced taste violations are relatively small, however, 
one could expand the fit function in powers of OLem = e 1 /{An). Thus, inclusion of a EM analytic 
terms to the fit function should allow the higher-charge data to be fit. That approach seems to work, 
and will be explored more in the future. For more details on EM taste- violating effects, see Ref. [Q]. 

Results and Outlook. - Figure ^| shows a typical fit of our data for AM 2 with physical quark 
charges to Eq. (||) (with added analytic NNLO terms). We fit partially quenched charged- and 
neutral-meson data simultaneously, but only the (unitary or approximately unitary) charged-meson 
data is shown in the plot. This fit has 55 data points and 26 parameters; other fits have as many 
as 120 data points, and from 20 to 30 fit parameters, depending on how many of the NNLO terms 
are included, and whether small variations with a 2 of the LO and NLO low-energy constants are 
allowed. The covariance matrix of the data is nearly singular, and the statistics are insufficient to 
determine it with enough precision to yield good correlated fits, so almost all fits currently used 
are uncorrelated. The fit shown has an (uncorrected) p value of 0.09. We note that what appear to 
be big discretization effects are actually due in large part to mistunings of the strange-quark mass, 
which is off by about 50% on the a = 0.12 fm ensembles and 25% on the a = 0.09 fm ensembles, 
but only by 2% on the 0.06 fm ensemble. 

The black and brown lines in Fig. || show the fit after setting valence and sea masses equal, 
adjusting m s to its physical value, and extrapolating to the continuum. The black lines adjust the 
sea charges to their physical values using NLO ZPT, while the brown line keeps the sea quarks 
uncharged. In the pion case, the adjustment vanishes identically, so no brown line is visible. In 
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Figure 1: A sampling of our partially quenched data in r\ units for EM splittings of pseudoscalar mesons 
with charge ±e p h ys , plotted versus the sum of the valence-quark masses. For clarity, only about a quarter of 
the data is shown. The red squares and magenta crosses show results for the two ensembles that differ only 
by the spatial volume: 20 3 and 28 3 , respectively. The vertical black bar labeled "BMW" shows the expected 
difference for kaons between these two volumes, based on the results from the BMW collaboration Next 
to it, the two points encircled in black are our "kaon-like" points for the volumes. 
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Figure 2: Typical SZPT fit to the squared-mass EM splitting AM 2 vs. the sum of the valence-quark masses. 
Only a small subset of the charged-meson data is shown. The red, blue and green curves correspond to three 
different lattice spacings. The brown and black curves are continuum limits for AM 2 , without or with the 
correction from XPT for physical sea-quark charges. The purple curves are the continuum limits for the 
K + -K° splitting (right), and the 7T + -"7r " splitting (left). 



the kaon case, the adjustment is a very small correction. From the black lines for the K + and 
K + , we subtract the corresponding results for the neutral mesons, "7T°" and K°, giving the purple 
lines. Results for (M 2 + — M 2 n0 „) r and (Mf+ - M^y are then obtained from the intersections of 
the purple lines and the vertical dashed-dotted lines that give the location of the physical point for 
each meson. The excellent agreement of the result for (M 2 + — M 2 ^,,) 7 and the experimental pion 
splitting (horizontal dotted line) is accidental, since our result has roughly 20% total error. 
We find the following preliminary results: 

(M|+ -M^f = 2100(90)(250) MeV 2 , (M 2 ) r = 901(8)(9)(?) MeV 2 

(M 2 + -M 2 ^,y = 1270(90)(230)(80) MeV 2 , (M 2 ^-) 7 = 157.8(1.4)(1.7)(?) MeV 2 

e = 0.65(7)(14)(10) (4) 

The first two errors in each case are statistical and lattice systematic uncertainties. The latter error 
comes largely from the effects of changing the assumptions entering into the chiral/continuum fit. 
Note, however, that finite- volume errors are not included at present. We expect that ultimately they 
will be a significant, but subdominant, source of error. The "?" for (M^y and (M 2 nQ „y represent 
the effect of EM quenching and, for (M 2 ^,,) 7 , the effect of neglected disconnected diagrams. These 
errors are likely to be much larger than the small quoted errors. For (M 2 + — M 2 ^,,) 7 and e the third 
error is a rough guess of the effect of neglecting disconnected diagrams, which we estimate by 50% 
of the result for {M 2 nl) „y. If we redefine e by replacing our computation of the pion EM splitting in 
Eq. ([]]) with the experimental splitting, we get e = 0.66(7) (20), which has larger chiral/continuum- 
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extrapolation errors, but no error from neglecting disconnected contributions. Our result for e is 
compatible with results from other groups |J. Using our values for s and (M^. ) r (assuming 
100% EM-quenching error in the latter quantity), our preliminary estimate for the EM uncertainty 
in m u /md is reduced by approximately a factor of two Jl0| ] from our previous error 

We are currently finishing the analysis of two more-chiral ensembles at a 0.06 fm in order 
to improve the chiral and continuum extrapolations, and are increasing the statistics for our finite- 
volume study. A "second-generation" calculation on HISQ ensembles, including ones at physical 
quark masses, is now beginning, and promises significant reductions in systematic errors. Calcula- 
tions with dynamical EM effects included are also being contemplated. 
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